Problemas Tema 4 


Ejercicio 4.1.- Sabiendo que las leyes fenomenológicas son lineales, deducir las unidades en 
el S.l. de la constante de proporcionalidad L de la ecuación (4.1) para cada uno de los 
procesos indicados en la Tabla 4.1, es decir, si la variable Y es (i) temperatura, (ii) cantidad 
de movimiento, (iii) concentración y (iv) potencial eléctrico. 

Las unidades de L, vendrán dadas, teniendo en cuenta (4.11) por Il| = -píL 


Propiedad | Propiedad | Variable Y 


| d Y/dx | 


energía 

J 

temperatura 

K 

K nr 1 

K tPr 

J K -1 nr 1 s ' 1 

impulso 

kg m s 1 

velocidad 

m s 1 

S' 1 

n 

kg nr 1 s ' 1 

materia 

moles 

concentración 

mol nr 3 

mol rrr 4 

D 

m 2 s 1 

carga 

culombio 

dif. potencial 

voltio 

Voltio m 1 

^rond 

rrr 1 Q 1 









Ejercicio 4.2.- Una celda cúbica de 0,100 m de lado se rellena con benceno. La cara superior 
se mantiene a 25 g C y la opuesta inferior a 15 g C. Calcular la cantidad de calor que fluye a 
través del benceno en una hora, una vez se haya alcanzado el régimen estacionario (sin 
convección). 


Para calcular la cantidad de calor necesitaremos conocer la conductividad térmica del 
benceno y el gradiente de temperaturas en la dirección del flujo que será la vertical dT/dz, 
pues el calor fluirá de la cara superior caliente a la inferior más fría. En la Tabla 4.2, vemos 
que la conductividad térmica para el benceno a 1 atm de presión y 22,5 Q C es k=0,1582 J K 1 
m^s' 1 , valor que a falta de mayor información consideraremos constante en el intervalo del 
temperaturas del ejercicio. Si la temperatura depende sólo de la vertical z, en régimen 
estacionario podemos hacer: 


dT _ AT 
dz Az 


1QK 

0,100 m 


= 100 Km -1 


Aplicando (4.12), la ley de Fourier monodimensional: 

= ^Q = _ Ak — 

dt At Az 

AQ = -Ak — At = -0,01 m 2 x 0,1582 J K 1 m 1 s 1 x 100 K m 1 x 3600 s =-569,5 J 
A z pues A=0,01 m 2 , At=lh=3600s. 









Ejercicio 4.3.- El coeficiente de viscosidad del agua líquida a 20 °C es 0,001002 kg-m-l-s 1 . En 
una conducción semejante a la mostrada en la figura 4.4 calcular la fuerza por unidad de área 
requerida para mantener la placa superior moviéndose a 0,250 m s 1 si la conducción tiene 
una profundidad de 0,500 m. 


z 



M 


0,5 m 


v 


La componente del gradiente de velocidad (en régimen estacionario) tiene un valor medio de 

—i = °’ 250ms = o,5005 -1 
dz 0,500 m 


Por lo tanto, haciendo uso de la ley de Newton monodimensional, (ecuación 4.14), 

4 dv x 
F_=-r¡A —- 

dz 

si nos piden la fuerza por unidad de área (será una presión, P ), o sea que: 


(P) = ^ = -n^ = {-(0,001002£gm~V)(0,500 s~‘)l = 5,01x10 ^kgm-'s" 2 = 
A dz 1 1 

= 5,01xl0~ 4 lVm~ 2 = 5,01xl0 _4 Pa 
















Ejercicio 4.4.- El agua fluye a través de un tubo de 42 cm de longitud y 5,20 mm de radio. 

Si la diferencia de presión entre dos puntos es de 0,050 atm y la temperatura es de 20 °C, 
determinar el volumen de agua que fluye cada hora. 

La ley de Poiseuille (ecuación 4.22), en su forma diferencial, nos da la relación entre el caudal 
de un fluido de viscosidad q que circula por una conducción cilindrica (cuyas dimensiones se 
especifiquen), y el gradiente de presión que lo impulsa a avanzar. Si el fluido es no 
compresible (el agua líquida en este caso), podemos usar la ecuación (4.23): 

dV Ttr 4 dP ^ dV m 4 AP 

o =— =-^ o = — =- 

dt 8 r) (-dz) dt 8 q (-Az) 

Az = 42 cm = 42xl0' 2 m 
r = 5,20 mm = 5,20xl0" 3 m 
q = 0,01002 kg m 1 s 1 
AP = 0,050 atm (101325 N m 2 atm' 1 ) 

Sustituyendo valores en la ecuación obtenemos: 

jt(0,00520 m) 4 (0,050atmxl01325Nm _2 atm _1 ) 3 

AV =--— -At = 12,44 m 

8(0,001002 kgm _1 s _1 )(0,420m) 







Ejercicio 4.5.- Cuando se establece el régimen estacionario en un flujo de materia a través 
de una superficie de 0,45 m 2 se observa que la cantidad de sustancia que fluye por minuto 
es de 5,65 moles de la misma. Si el gradiente de concentración de dicha sustancia es 
7,25xl0' 2 M nr 1 , determinar el coeficiente de difusión de la sustancia en dicho medio 
disolvente. 

Para calcular el coeficiente de difusión, D, haremos uso de la Primera Ley de Fick 
(ecuación 4.28) que nos proporciona el flujo de materia a través de una superficie A, si se 
establece un gradiente de concentración en régimen estacionario: 


dn 

di 


= -D A 


de 

dz 


D = 


A n 1 

A t A(dc / dz) 


5,65 mol 


1 


lminx60smin 1 


0,45ra 2 x7,25xl0' 


mol 


3t-1 


LxlO m L 


m 


2,886-10 ~ 3 m 2 s~ l 







Ejercicio 4.6.- En textos de electricidad y electromagnetismo, aparece la ley de Ohm como 
la siguiente relación: V=IR. En los fenómenos de transporte hemos visto, Tabla 4.1 y 
ecuación (4.34), que la ley de Ohm se definía como J=-k V<£. 

Definir cada uno de los símbolos que aparecen en ambas leyes de Ohm y encontrar la 
relación entre los parámetros k y R. 


(a) Definición de los símbolos que se explicitan en el enunciado: 

- Ley de Ohm de textos de electromagnetismo: V= IR 

(1) V es el voltaje o diferencia de potencial entre dos puntos de un hilo conductor por el que pasa una corriente 
eléctrica. La unidad de medida es el voltio. 

(2) I es la intensidad de corriente eléctrica que pasa a través de la sección del conductor y se refiere a la cantidad de 
carga eléctrica que por unidad de tiempo atraviesa la sección del conductor. La unidad de medida es el amperio. 

(3) R es la resistencia que ofrece el conductor al paso de dicha intensidad. La unidad de medida es el ohmio. 

- Ley de Ohm de fenómenos de transporte: J=-kVc|d 

(1) J es la densidad de flujo de carga o densidad de corriente eléctrica y determina la cantidad de carga que circula 
por el conductor por unidad de área y tiempo por lo que se le asimila con una intensidad eléctrica por unidad de 
área. La unidad de medida es el amperio dividido por metro cuadrado. 

(2) k es un coeficiente de proporcionalidad del transporte de carga. Relaciona el transporte de la carga entre dos 
puntos en los que se ha establecido un gradiente de potencial. Se le denomina conductividad eléctrica o 
conductancia específica. La unidad de medida es el Siemens dividido por metro. 

(3) Vcf) es el gradiente de potencial eléctrico y mide la diferencia de potencial eléctrico entre dos puntos por 
unidad de separación en el espacio de los mismos. La unidad de medida es el voltio dividido por metro. 


(b) Relación entre los parámetros k y R 


V diferencia de potencial 
I int en si dad de corriente 

J densidad de corriente eléctrica 
V(|) gradiente de potencial 

intensidad distancia entre 
diferencia de potencial sección 


int ensidad 

= _ sección _ 

diferencia de potencial 

distancia entre puntos 
puntos 


Si llamamos Z 

a la distancia entre puntos (o sea la longitud del conductor) y S a la sección del mismo, es obvio que comparando 
ambas expresiones anteriores se tiene que: 

_ í i _ i l _ l 

k "vs"rs"rs 

Entre la conductividad eléctrica y la resistencia existe una analogía recíproca tal que se suele definir un nuevo 
parámetro llamado resistividad específica, p, que depende de la naturaleza del conductor, la temperatura y la presión, 
definido como: 

s R 

con lo que 


s 

i i 

a menudo nos encontramos con la relación recíproca entre k y p, es decir: K = — 

P 



Ejercicio 4.7.- Despreciando la diferencia de masas, calcular el coeficiente de difusión de 
moléculas de nitrógeno datado isotópicamente en nitrógeno ordinario a 298K y 1 atm de 
presión. Dato: diámetro molecular del nitrógeno ordinario, d = 3,7x10 10 m. 

El coeficiente de difusión, D^*, en concreto viene relacionado con los parámetros de la TCG 
por la ecuación 4.68: 



Sustituyendo los valores conocidos en la ecuación 4.68 se tiene: 


3 1 (8,3145 JK^mor 1 ) 3 * (298 K) 3 _1_ 

jt(0,028kgmol“ 1) (6,02214xl0 23 mor'HSJxKr'V) 2 101325Nm" 2 


I 


= l,87xl0' 5 m 2 s _1 












Ejercicio 4.8.- La conductividad eléctrica, k del agua pura es 5,5xl0' 6 Q 1 m 1 a 25°C. ¿Cuál 
es el valor del producto iónico del agua, K w =[H + ][OH]? 

La ecuación (4.48) nos relaciona la conductividad eléctrica, k, con la concentración 

k = Fc ( V+ u 0H _) 

por lo que podemos obtener la concentración, c, de uno de ellos (que será igual a la del 
otro) a través de: 


K 


5,5xlO^Q —1 m~ 1 


F(u H+ +u 0fJ _) 96485Cmol" 1 


) 96485Cmor 1 (36,25 + 20,64)xl(T 8 m 2 V~ 1 s' 1 


= l,002xl0~ 4 ^ = l,002xl0 _7 molL _1 


m 







Ejercicio 4.9.-Dos depósitos de calor con temperaturas respectivas de 325 y 275 K se ponen 
en contacto mediante una varilla de hierro de 200 cm de longitud y 24 cm2 de sección 
transversal. Calcular el flujo de calor entre los depósitos cuando el sistema alcanza su estado 
estacionario. La conductividad térmica del hierro a 25 °C es 0.804 J K-l cm-1 s-1.(Solución: 
4.824 J s-1) 



, 1 dQ dT 

A dt dz 


dQ A dT 

= -kA- 


dt 


dz 


Para calcular la fuerza (gradiente de T con z) podemos utilizar el hecho de que al 
alcanzar el estado estacionario tendremos un perfil lineal: 



dT 

dz 


AT 

Az 


T 2 -Ti 


Z 2 -Zi 


275-325 

200 


= -0.25K-cm 


-1 


Sustituyendo en la ley de Fourier 


— = -k-A-— = -(o.804JK- 1 cm- 1 -s- 1 }(24cm 2 )(-0.25K-cnr 1 )= 4.824J-S' 1 

Aunque hemos mezclado unidades, nótese que todos los cm se van, quedando 
unidades del SI. 

El resultado del flujo es positivo, lo que indica que el calor va del foco caliente al 
frío, como debe de ser 



























Ejercicio 4.10.- Calcular la conductividad térmica del He a 1 atm y 0 °C y a 10 atm y 100 °C. 
Utilizar el valor del diámetro molecular que se obtiene a partir de medidas de viscosidad a 1 
atm y 0 °C, d = 2.2 Á. El valor experimental a 0 °C y 1 atm es 1.4-10-3 J K-l cm-1 s-1.(Solución: 
1.42110-3 y 1.66 10-3 J cm-1 K-ls-1) 




P = 


X = 


v = 


K = 


25 / RT \ 


1/2 


32 y jtM 


N A d 2 ° v ’ m 


N _ P 

V “ kT 
1 kT 

V2jtd^ P 

1/2 

/ 8kT \ 

\ jtm / 


d=2.2-10 10 m 
M=4.003-10' 3 kg-mol 1 
C v m =3/2R (gas monoatómico) 


Datos que tenemos: 





















a) T=273.15 K y P= 1 atm=101325 Pa 


25 / RT \ 1/2 

K ” 32 1, itM j N A d 2 ° vm 

25/ 8.31451JmolK _1 -273.15K\ 12 1 [3.,,.,., 

32 y Jt4.003-10“ kg-mol“ ) 6.02214-10“ 23 mol“ 1 (2.2-10“ 1 °m)^ [2 

= 0.142Jm" 1 K _1 s' 1 = 1.42-10 _3 Jcm' 1 K _1 s' 1 
b) T=373.15 K y P= 10 atm=1013250 Pa 


25 / RT \ 
32 ^) 


-^C„ m =0.166 J-m _1 -K _1 -s' 1 = 1.66-10' 3 J-crrf 1 K _1 -s _1 
N A d 2 vm 


1/2 



Ejercicio 4.11.- La viscosidad y la densidad de la sangre humana a la temperatura del cuerpo 
son 4 cP y 1.0 g cm-3, respectivamente. El flujo de la sangre desde el corazón a través de la aorta 
es 5 L min-1 en un cuerpo humano en reposo. El diámetro de la aorta es típicamente de 2.5 cm. 
Calcule: (a) el gradiente de presión a lo largo de la aorta; (b) la velocidad media de la sangre; 


Ecuación de Poiseuille para líquidos 

r]= 4 cp= 4-10 3 N s nrr 2 
r= 1.25-10- 2 m 

AV/At= 5 l-min'^ (5/60)-10" 3 nrP-S" 1 


AV 

jtr 4 

AP 

At " 

8r| 

l 


AP 8r| AV 

i jtr 4 At 


-34.77N-rrf 3 


-34.77Pa-m' 1 


b) Velocidad media 

En un tiempo t el fluido avanza (en promedio) una distancia <d>, fluyendo un volumen V 


<d> 



Area S 













Entonces, la velocidad de flujo o volumen que circula por unidad de tiempo será 


V _ < d >-S 

T“ t 

La distancia media recorrida por unidad de tiempo es la velocidad media: 


Quedando: 


V/t 

< V >=- 

S 



(5/60)-10' 3 m 3 s' 1 

jt(1.25-10' 2 m) 2 


0.17m-s _1 



Ejercicio 4.12.-Dos tubos de cobre, cada uno de 3 m de longitud, con un diámetro interno el 
primero de 2.6 cm y de 1.3 cm el segundo, se conectan en serie. Se establece una presión de 5 
atm en el extremo abierto del tubo más ancho, y del extremo más estrecho sale aceite a una 
presión de 1 atm. Para el aceite, h = 0.114 Pa s a 15 °C. a) Calcule la presión en el punto en que 
se unen los dos tubos, b) ¿Cuántos litros por minuto pueden obtenerse mediante esta 
combinación? 



P f =1 atm 


Cuando se conectan 2 tuberías se cumple que el volumen que circula por unidad de tiempo 
es igual en ambas A . . á ^ ^ 

AV jir, P u - P| 

Para el tramo 1 podemos escribir: — = - 

At 8r| z u - z¡ 

AV nr_ 4 P f - P, 

Para el tramo 2 podemos escribir: - 


At 


8r| z f - z u 


Igualando y teniendo en cuenta que los dos tramos miden kjual: 


jtr, 


P -P 

U I 


< P f - P u _ 

8ri z f - z u 8r] z u - z¡ 



V r i / 


P -P 

U I 

p -P 

r f ' u 





















La única incógnita que queda es la presión en el punto de unión P u 


/IV P.-5 


v2/ 


1-P. 


* 


P=4.765 atm 


El perfil de presiones que tenemos es por tanto: 


Tramo 1 


\Tramo 2 


u 


El flujo se iguala porque la caída de presión se reparte de manera desigual entre los 2 
tramos. El flujo lo podemos calcular usando cualquiera de los tramos: 


AV JCTo Pf — P.j JlT P — P¡ QO , a _4 3-1 a n noi 

— = — 2 f ^ = — 1 U L = 7.82-10 rrrs = 46.921-min 

At 8r| z f - z u 8r| z u - z¡ 













Ejercicio 4.13.- La viscosidad del 02 a 0 °C y presiones del orden de magnitud de 1 atm es 
1.92-10-4 P. Calcular el flujo, en g-s' 1 , del 02 a 0 °C a través de un tubo de 0.420 mm de 
diámetro interior y 220 cm de longitud, cuando las presiones a la entrada y salida son de 2.00 y 
1.00 atm, respectivamente. 

La ecuación de Poiseuille en forma diferencial es: _ = 

dt 8r| dz 


Para gases no se puede integrar directamente, puesto que el volumen es función de la presión. 
Para integrarla podemos expresar el flujo en masa que circula por unidad de tiempo: 

dV RT dm 
~dF " PM dt 

Sustituyendo en la ec. de Poiseuille nos queda: 

RT dm jtr 4 dP 

PM dt 8r| dz 

En el régimen estacionario, la masa de gas que circula por unidad de tiempo es una constante 
por lo que podemos integrar la ecuación anterior: 


Am , jrr 4 M 

— dz =-PdP 

At 8r|RT 



Am 

jtr 4 M P f 2 -P¡ 2 

At ~ 

16r]RT z f -z¡ 



Ecuación de Poiseuille para 
gases 



















Sustituyendo los datos que nos da el problema: 


Am jir 4 M P 2 - P 2 


At 16r|RT z f -z¡ 

jt(2.1-10 _4 m) 4 (32-10' 3 Kg-mol' 1 ) (I 2 


16(1.94-10“ 5 Pas)(8.3145J-K' 1 mol' 1 )273.15K 


— = 3.88-10' 6 Kgs _1 = S^S-IO^g-s" 1 
At 


2 2 )-(l01325Pa) 2 

2.20m 



Ejercicio 4.14.- Calcule la velocidad final de caída de una bola de acero de 1.00 mm de 
diámetro y 4 mg de masa, en agua a 25 °C. Repita el cálculo para glicerina (densidad 1.25 g 
cm-3). Las viscosidades del agua y de la glicerina a 25 °C y 1 atm son 0.89 y 954 cP. 
respectivamente. 

Cuando la bola cae en el interior de un fluido hay tres fuerzas actuando sobre ella: el peso, el 
empuje y el rozamiento (relacionado con la viscosidad del medio y la velocidad de la bola) 



v P 

Si el cuerpo es más denso que el fluido, el peso es mayor que el empuje y la bola cae en su 
interior. Al haber una fuerza resultante, la bola se irá acelerando poco a poco. Sin embargo, a 
medida que aumenta su velocidad, aumenta también el rozamiento, con lo que llega un 
momento en el que la suma de todas las fuerzas se anula. Es la condición de estado 
estacionario, momento en el que la bola pasa a moverse con velocidad constante 


El peso, m g, se puede escribir teniendo en cuenta el volumen y densidad de la bola: 

n 4 3 Siendo p h la densidad de la bola 

P = mg = ,tr 3 p b g 

El empuje es el peso del volumen de fluido desplazado: 





Por último, la fuerza de rozamiento de un cuerpo esférico de radio r en el interior de un fluido de 
viscosidad r\ viene dado por la ley de Stokes, y es función de la velocidad v con la que se 
mueve el cuerpo en el fluido 


F t] = 6jtr|rv 


En el estado estacionario se cumplirá: 


2 F -°. 


P-E-F -0 


11 


Sustituyendo las expresiones anteriores: 

| ™- 3 p b g -1 Jtr 3 p f g - 6jtrr|V = 0 

Simplificando y despejando la incógnita (v): 

2gr 2 (p b -p f ) 


V = 


9r| 


Esta relación puede emplearse también para, una vez medida expenmentalmente la velocidad 
de caída, calcular la viscosidad del fluido. 



a)agua 


p b =m/V=7.64 g-crrr 
p f =1.0 g-crrr 3 
t]= 0.89-10' 2 Posies 
g=980.665 cms 2 


b) glicerina 

p b =m/V=7.64 g-crrr 
p f =1.25 g-crrr 3 
r]=9.54 Posies 
g=980.665 cm s 2 


2 sr 2 jPb-P f ) 

9rj 


406cm\s' 1 


2gr 2 (Pb Pf) . o.36cm-s _1 
9r| 






Ejercicio 4.15.- ¿Con qué velocidad, pueden ascender las burbujas de aire (cavidades) en 
agua a 25 °C si sus diámetros son de 1 mm? Datos adicionales del agua a 25 Q C: densidad, 10 3 
kg-nr 3 , viscosidad 8.91-10' 4 kg-rn^s 1 


En este caso las fuerzas actúan de acuerdo con el siguiente esquema: 



YP = 0, E-P-F 11 =0 


Y 


p 


Sustituyendo : 


|jtr 3 p f g -1 Jtr 3 p aire g - 6jirriv = 0 

Simplificando y despejando la incógnita (v): 


2gr 2 (p f - Paire) 2 9 r2 Pf 
9r| 9 t] 

2r 2 Pf g 2(5x10~ 4 m) 2 (1x10 3 kg m~ 3 )(9,806 m s~ 2 ) i 

9(8,91x 1 O -4 kg m _1 s _1 ) 


61,1cms 1 











Ejercicio 4.16.- Las viscosidades del C02(g) a 1 atm y 0, 490 y 850 °C son 139, 330 y 436 \iP, 
respectivamente. Calcule el diámetro de esfera rígida aparente del C02 a cada una de estas 
temperaturas. 


5 (MRT) 172 

^ 16 Vjt N A d 2 

Despejando el diámetro: d 2 


5 (MRT) i;2 
16 n N A r| 


Los datos son: M=44-10‘ 3 kg mol‘ 1 

R=8.3145 J K'^mol' 1 
N A =6.022-10 23 mol' 1 


a) T=273.15 K y r|=139 10 7 Pas 

b) T=763.15 K y ri=330-10" 7 Pa s 

c) T= 1123.15 K y ri=436-10' 7 Pa s 


d 2 =2.106'10" 19 m 2 „ d= 4.59-10- 10 m = 4.59 A 
d= 3.85-10" 10 m = 3.85 A 


d= 3.69-10" 10 m = 3.69 A 


En principio, el diámetro de esfera rígida debería ser constante, con lo que estos resultados 
nos muestran las limitaciones de esta aproximación. 

El diámetro disminuye porque al aumentar la temperatura aumenta la velocidad de las 
moléculas, pudiendo producirse un mayor acercamiento de las mismas durante la colisión. 








Ejercicio 4.17.- El hidrógeno gaseoso se difunde a través de una lámina de paladio de 0.0050 cm de espesor. Del lado 
izquierdo de la lámina, el hidrógeno se mantiene a 25.0 °C y una presión de 750 mm, mientras que del lado derecho se 
mantiene un buen vacío. Después de 24 h, el volumen de hidrógeno en el compartimento de la izquierda disminuye en 
14.1 cm3. Si el área de la lámina a través de la cual ocurre la difusión es 0.743 cm2. ¿Cuál es el coeficiente de difusión del 
hidrógeno en el paladio? 



5-10 -3 cm 



Si mantenemos las presiones de hidrógeno constantes 
se alcanzará un estado estacionario con un perfil lineal 
Paladio D 

*1 /m o a - 1_3 


Ci = 


c 2 - 


i T 

RT 


= 40.34mol-m 


= 0 


a cada lado de la lámina entonces 
de concentraciones en la lámina de 


Con lo que el gradiente de concentraciones será 


de 

dz 


Ac 0-40.34 
Az ~~ 5-10' 5 


-8.067-10 5 molrrf 4 






















El coeficiente de difusión lo podremos obtener de la primera ley de Fick: 


dn 

= -D A dC 

dt 

dz 


Para despejar el coeficiente de difusión necesitamos saber lo que vale el flujo, que al 
alcanzarse el estado estacionario lo podemos obtener simplemente como el número de 
moles que han pasado de una lado al otro dividido por el tiempo. El número de moles que 
han pasado de la izquierda a la derecha lo podemos obtener por la disminución de volumen 
que se ha producido, ya que la presión en ese lado permanece constante: 


PAV (750/760)101325-14.1-10' 6 

dn = An = = 8.3145298.15 

dt At At 24-3600 

Quedando para el coeficiente de difusión 


6.583-10 _9 mol-s _1 


dn 

D = dt 

A — 
dz 


-6.583-10' 9 mol-s _1 
0.743-10" 4 m 2 ■ (-8.067-10 5 m ol -m -4 ) 


= 1.098-10 _1 ° m 2 -s _1 















Ejercicio 4.18.- El diámetro molecular que se obtienen para el 02 a partir de medidas de 
viscosidad a O °C y 1 atm es 3.6 Á. Calcular el coeficiente de autodifusión del 0 2 a O °C y 
presiones de 1.00 atm y 10.0 atm. El valor experimental a 0 °C y 1 atm es 0.19 cm 2 s' 1 

Para calcular el coeficiente de autodifusión usamos la expresión proporcionada por la teoría 
cinética de gases (versión rigurosa) 




3 íkT \ í/2 kT 

D - 'k(v) 

16 


8d 2 \ nm ) P 


a) Datos: T=273.15 

m=32 -10' 3 /N A kg 
P= 1 atm= 101325 Pa 
d=3.6-10' 10 m 


D=1.62 ■ 10-Ws- 1 =0.162 cm 2 S ' 1 
El error es del 15% aprox. 


a) Datos: T=273.15 

m=32 -10 _3 /N A kg 
P= 10 atm= 1013250 Pa 
d=3.6-10' 10 m 


D=1.62-10- 6 m 2 s- 1 =0.0162 cm 2 ^ 1 













Ejercicio 4.19. -Suponga un sistema unidimensional que se extiende desde z = 0 a z = oo. En el instante t = 0 hay N 0 partículas 
en el punto z = 0. Supuesta válida la segunda ley de Fick se ha deducido que: 


c(z,t) = 




(jrDt) 


- 2 — e 4Dt 


n 


Calcule cuál es la probabilidad de encontrar una partícula en una posición comprendida entre z y z+dz. Por último, 
calcule los valores de <z> y <z 2 >. NOTA:La concentración en un sistema unidimensional viene dada en "partículas por 
unidad de longitud". 

N 0 


0 

Probabilidad de encontrar una molécula entre z y z+dz en instante t será: 


dp(z,t) = 


dN(z,t) 

N„ 


El número de moléculas dN(z,t) se puede calcular como concentración por longitud. Teniendo en 
cuenta que es un sistema unidimensional: 


dN(z,t) = C(z,t)dz = 


K 


(jcDt) 1 ' 2 


Z 

4Dt dz 


Con lo que la probabilidad de encontrar una molécula entre z y z+dz en el instante t será: 

dN(z,t) C(z,t)dz 1 


dp(z,t)« 


N, 


N, 


(jtDt) 1/2 


z 

4Dt -dz 








Para calcular cualquier propiedad promedio hacemos uso de la probabilidad. Así, para <z>: 


OO 00 A Z A 00 Z 

<z>-f z dp( z ..).rz„e-=dz. j - 5r rze-=.dz 


Los límites (0,o°) vienen dados por el sistema que estamos estudiando. La integral se resuelve con 
ayuda de las tablas: 


< z >= 


(jtDt) 


1 00 1 1 
_ f 7'g 4Dt -Hv = _■_ 

»t) 1/2 J («Dt) 1,2 2 M 

4Dt 


De igual modo podemos operar para calcular <z 2 >: 


<z 2 >= Jz 2 dp(z,t) = 


(jrDt) 1 


'z 2 -e 4Dt -dz = 


(jtDt) 1 


= 2Dt 


1 

4Dt 



Ejercicio 4.20.-Una disolución concentrada de 10 g de sacarosa en 5 mL de agua se introdujo en un cilindro de 5 cm de 
diámetro. Posteriormente, se añadió un litro de agua con sumo cuidado para no perturbar la superficie de la capa de 
disolución. Calcule la concentración a 5 cm por encima de la capa transcurrido un tiempo de (a) 10 s y (b) 1 año. Ignore los 
efectos gravitacionales y considere únicamente el proceso de difusión. El coeficiente de difusión de la sacarosa a 25 °C es 
5.2 -10' 6 cm 2 s 1 . La solución de la 2a ley de Fick para este 2 caso es: 

c(z,t) = 


ADt 



A{jrDtf 2 

La concentración de sacarosa en función de z y t viene dada por: 


; (z,t) = 


n 


o 


_z 


4Dt 


z=0 


A(jrDt) 


(i) 


n 0 es el número de moles, que podemos calcular sabiendo la masa molar de la sacarosa 
(¿ 12 H 22 O 11 ). Si expresamos todos los datos en el sistema CGS: 


N 0 = 10 / 342.3 = 2.92-10 2 moles 
A=jcr 2 = 19.64 cm 2 
D= 5.2-10 6 cm 2 -s 1 

La concentración (en moles/cm 3 ) en z=5cm vendrá dada en función del tiempo sustituyendo los 
datos en (1): 





n 


o 


_z 


4Dt 


2.92-10 


-2 


25 


A(jrDt) 


4-5.2-10" 6 t 


-6 i .\ 1/2 


19.64-(jt-5.2-10t) 


0.368 

^1/2 


1 . 202 - 10 6 

e ' (mol/cm 3 ) 


(a) t = lOs 



0.368 

j.1/2 


1 . 202 - 10 6 

e 1 


= 0.1164-e 


- 1 . 202 - 10 5 



(b) t = 1año=3.1536-10 7 s 


c(5,3.1536-10 7 ) 


0.368 

^1/2 e 


1 . 202 - 10 6 

t~~ 


= 6.554-10' 5 -e 


= 6.3-10“ 5 ^ = 0.063M 
cm 


0.038 












Ejercicio 4.21.- Calcular la distancia cuadrática media recorrida por una molécula de glucosa en agua a 25 °C en 30 minutos. 
Suponer que las moléculas de glucosa se difunden a partir de (a) una capa depositada en el fondo del vaso y (b) un pequeño 
terrón suspendido en el seno del agua. ¿Cuánto tiempo tardarán las moléculas de glucosa en recorrer una distancia de 1 
mm y 1 cm desde su punto de partida en el caso a? El coeficiente de difusión de la glucosa en agua a 25 °C es 0.673-10-9 m2 
s-1. Las soluciones de la 2a ley de Fick son: 


a) c(z,t) = 


Un 


__Z _ 

4 Dt 


A(jrDt) 


¡12 ' 



b ) c(z,t) = 


n, 


r~ 

ADt 


8 (jiDt) 


3¡2 


oo oo 

<z 2 >= Jz 2 -dp(z,t)= Jz 2 -f(z,t) dz 


o o 

La probabilidad de encontrar un mol 


Siendo f(z,t) la función de distribución 
de azúcar entre z y z+dz en el instante t será: 


dp(z,t) = 


dn(z,t) _ C(z,t) A dz _ 1 

n o n o (jtDt)' 


Z 

401 dz 


Por otro lado la probabilidad se puede escribir como: 


dp(z,t) = f(z,t)dz 











Así, el valor medio de z 2 será: 



Y por tanto 

< z 2 >= 2Dt = 20.673-10' 9 m 2 s' 1 -1,8-10 3 s = 2.423-10“ 6 m 2 

Y la raíz de la distancia cuadrática media: 

z rms = (< z2 >) 1,2 = 1.56-10" 3 m 


¿Cuánto tiempo tardarán las moléculas de glucosa en recorrer una distancia de 1 mm y 1 cm desde su 
punto de partida en el caso a? 


2 2 

j- _ < ^ ^ _ ^rms 


z =10 3 m 

rms 


z = 10 2 m 

rms 


t = 743s 


2D 2D 


t = 74294s = 20.6h 











00 



< r 2 >= Jr 2 f(r,t)dr 

0 


Probabilidad de encontrar un mol entre r y r+dr en instante t: 


dp(r,t) = 


dn(r,t) 


n, 


Quedando la función de 
distribución como: 


n, 


n, 


n, 


rinfr n dn(r,t) c(r,t)dV c(r,t)4itr 2 dr 4jtr 2 ,y 
P( ’ ,= = = = 8(^Dt) 3 ^ 

dp(r,t) = f(r,t)dr 

1 



< r 2 >=J r 2 f(r)dr = Jr 

O O 


2~- r / 4 Dt 


re 


2jr 2 (Dt) 


dr = 


2jt /2 (Dt) : 


■p 

o 


r 4 e ,/4Dt dr = 6Dt 


Y por tanto 

< r 2 >= 7.268-10~ 6 m 2 

Y la raíz de la distancia cuadrática media: 

W = (<r 2 >) 12 =2.70-10“ 3 m 

















Ejercicio 4.22.- El coeficiente de difusión del níquel en cobre es 10' 9 cm2 s-1 a 1025 °C. 
Calcular el tiempo necesario para que los átomos de níquel se difundan una distancia de 1 cm 
en el cobre. Repita el cálculo para la difusión del aluminio en cobre a 20 ° 


Si consideramos la difusión en 1-D 

< Z 2 >= 2Dt j- _ Z rms 

2D 



2 2 

Ni en Cu t = Zfms = - ^ = 5-10 8 s = 15.85años 

2D 2-10' 9 

2 f 

AL en Cu t = = _ n = 5-10 29 s = 1.6-10 22 años 

2D 2-10' 30 







Ejercicio 4.23.-Estimar el tiempo requerido por las moléculas de un neurotransmisor para 
difundirse a través de una sinapsis (separación entre dos células nerviosas) de 50 nm, si su 
coeficiente de difusión a la temperatura del cuerpo humano es 5 x 10-10 m 2 s -1 . 



Si consideramos la difusión en 1-D 

< Z 2 >= 2Dt [ — Z rms 

2D 




2D 


(5-1 o- 8 y 

2-5-1 (T 10 


2.5-10“ 6 s 







Ejercicio 4.24.- La gutagamba es una resina gomosa que se extrae de árboles originarios de 
la selva de Camboya. Las observaciones de Perrin sobre partículas esféricas de esta resina, 
de un radio medio de 2,1 -10' 5 cm en suspensión acuosa, a 17 °C ( r¡ - 0,011 P), condujeron 
a los siguientes resultados para los valores de z rms : 7,1 -10' 4 , 10,6 -10' 4 y 11,3-10' 4 cm, para 
intervalos de tiempo de 30, 60 y 90 s , respectivamente. A partir de estos datos, calcular el 
número de Avogadro. 


Cuando se estudia la difusión de partículas esféricas en medios viscosos, el 
desplazamiento cuadrático medio, para un tiempo, t, puede expresarse por la fórmula de 
Einstein (ecuación 4.102): 

z 2 = k ¿ r t 
3jtr|T 


En dicha expresión, todas las magnitudes son medióles, excepto k B , que puede deducirse 
de la misma haciendo: 



3jrrf z 2 

Tt 


3jtr|rz 


2 

rms 


Tt 


Finalmente y como en la época de Perrin ya se conocía el valor de R, (la constante de los 
gases), 


N 


A = 



RTt 

3jtr|rz rms 


= o 


t 

^-rms 







El producto de constantes, será: 


o = 


RT 

3jrr|r 


8,31451 JK~ 1 mol~ 1 x290,15K 
3 it x 1,1x10 -3 Pas x 2,1x10 -7 m 


1,108x10 12 m 2 s _1 mor 1 


a) z rms = 7,1-10-4 cm = 7,1-10-6 m ; t= 30 s 

El valor del número de Avogadro, será: 

N A (a) = 1,108x10 12 m 2 s“ 1 mor 1 7 -——^ = 6,59x10 23 mol' 1 

7,1x10“ 6 mf 


b) z rms = 10,6-10-4 cm = 10,6-10-6 m ; t= 60 s 

El valor del número de Avogadro, será: 

N a (b) = 1,108x10 12 m 2 s“ 1 mor 1 7 -——^ = 5,91x10 23 mol' 1 

(i 0,6x10“ 6 mf 

c) z rms = 11,3-10-4 cm = 11,3-10-6 m ; t = 90s 

El valor del número de Avogadro, será: 

N a (c) = 1,108x10 12 mV 1 mor 1 7 -——^ = 7,81x10 23 mol' 1 

(i 1,3x1 0 -6 mjr 



Como se observa, la dispersión es muy grande. Las incertidumbres en la determinación de 
un radio medio pueden ser cruciales. No obstante el valor medio entre los tres resultados 
obtenidos es de 6,77'1023 que tiene un error del orden del 12 % respecto del valor 
aceptado actualmente (N A = 6,02214 10 23 mol 1 ). 



Ejercicio 4.25.- Calcular el coeficiente de difusión para una molécula de hemoglobina (masa 
molecular 63000, d = 50 Á) en disolución acuosa a 20 °C. Comparar su resultado con el valor 
experimental de 6,9x1o 11 m 2 s' 1 a 20 °C y con el valor que se obtendría si fuera válida para los 
líquidos, la aproximación de colisiones binarias de los gases (d=3.2 Á) 


Cuando se estudia la difusión de partículas esféricas en medios viscosos, el coeficiente de 
difusión se expresa mediante la ecuación de Stokes-Einstein: 


n 00 

u He,H 2 0 


k B T 

6 ™lH 2 O r He 


k B = 1,38066'1023 J K'i 
T =20 °C = 293,15 K 
r¡ - 0,01 P = 10-3 Pa-s 
r He = 2,5 Á = 2,5-10 10 m 


D 


00 


He,1-1,0 


1,38066x1 0~ 23 x 293,1 5 
6x jtx10 _3 x25x10 _10 


= 8,59x10 _1 Ws -1 


El error relativo respecto del valor experimental es: 

8,59 - 6,9 


e r = 


0,245 = 24,5% 


6,9 






(b) Supongamos que queremos calcular la viscosidad del agua a través de la fórmula derivada 
de la teoría cinética de los gases (ecuación 4.63): 

5npX v 5 (MRT) 1/2 

*' H ° 32 16 ji N A d^ 0 

Sustituyendo los datos, tenemos: 



5 (18x10~ 3 x8,31451x293,16) 1/2 

16 n 6,02205x10 23 x(3,2x10" 10 ) 2 


1,89x10 5 Pas 


Si calculamos ahora el coeficiente de difusión de la Hemoglobina usando esta viscosidad para 
el agua: 


D’= D°—=8,59x10 _11 x 


-3 




10 


1,89x10 


-5 


= 4,54x10“ 9 m 2 s _1 


Como puede verse este resultado es irreal y es consecuencia de la no validez de las 
hipótesis de la teoría cinética de los gases para los líquidos. 



Ejercicio 4.26.- La constante de difusión de la hemoglobina en agua a 20 °C / es 6,9-10 11 
m 2 s 1 y la viscosidad del agua a esa temperatura es 1,002-10' 3 kg-m^s -1 . Suponiendo que las 
moléculas son esféricas, calcular el volumen molar de la hemoglobina y compararlo con el 
valor experimental (la densidad de la hemoglobina a 20 °C es 1,335 g-mL 1 ). 


r He = 


k B T 

6 ™lH 2 0 D He,H 2 0 


Datos : k B = 1,38066-10' 23 J K 1 
NA = 6,022 -10 23 mol 1 
T =20 °C =293,15 K 
r¡ = 0,01 P = lxlO' 3 Pa s 
D = 6,9-10 11 m 2 s' 1 

M = 63000 g mol' 1 
r = 1,335 g mL 1 = 1,335-10 3 g L 1 


1,38066x10~ 23 x293,15 
6ji x 1x10 -3 x6,9x10 _11 


= 3,11x10 -9 m = 31,1Á 





El volumen molar experimental es fácil de obtener a partir de la densidad: 


-i 


V = 


Mgmol 
P gL 1,335x10" 3 


63000 , 7 , n , ,_i 

= 47,19Lmol 


Si consideramos las moléculas de hemoglobina como esferas, tendremos: 


Amolar - v 

con lo que ahora: 

Vmolar = |x nx (3,1 1 x 10 ‘ 9 f x 6 , 022 x 10 23 = 0,0758m 3 mor 1 = 75,8Lmol 


-1 


En este caso, el volumen teórico es muy superior al que se deduce de la densidad de la 
hemoglobina. Esto es una prueba de que la hemoglobina ocupa un volumen inferior al que 
correspondería como esfera para su radio aparente. La molécula de hemoglobina no es 
esférica. 





